Variable paralog expression underlies phenotype variation.

Human faces are variable; we look different from one another. Craniofacial disorders further increase facial variation. To understand craniofacial variation and how it can be buffered, we analyzed the zebrafish mef2ca mutant. When this transcription factor encoding gene is mutated, zebrafish develop dramatically variable craniofacial phenotypes. Years of selective breeding for low and high penetrance of mutant phenotypes produced strains that are either resilient or sensitive to the mef2ca mutation. Here, we compared gene expression between these strains, which revealed that selective breeding enriched for high and low mef2ca paralog expression in the low- and high-penetrance strains, respectively. We found that mef2ca paralog expression is variable in unselected wild-type zebrafish, motivating the hypothesis that heritable variation in paralog expression underlies mutant phenotype severity and variation. In support, mutagenizing the mef2ca paralogs, mef2aa, mef2b, mef2cb, and mef2d demonstrated modular buffering by paralogs. Specifically, some paralogs buffer severity while others buffer variability. We present a novel, mechanistic model for phenotypic variation where variable, vestigial paralog expression buffers development. These studies are a major step forward in understanding the mechanisms of facial variation, including how some genetically resilient individuals can overcome a deleterious mutation.

The alx3 gene shapes the zebrafish neurocranium by regulating frontonasal neural crest cell differentiation timing.

During craniofacial development, different populations of cartilage- and bone-forming cells develop in precise locations in the head. Most of these cells are derived from pluripotent cranial neural crest cells and differentiate with distinct developmental timing and cellular morphologies. The mechanisms that divide neural crest cells into discrete populations are not fully understood. Here, we use single-cell RNA sequencing to transcriptomically define different populations of cranial neural crest cells. We discovered that the gene family encoding the Alx transcription factors is enriched in the frontonasal population of neural crest cells. Genetic mutant analyses indicate that alx3 functions to regulate the distinct differentiation timing and cellular morphologies among frontonasal neural crest cell subpopulations. This study furthers our understanding of how genes controlling developmental timing shape craniofacial skeletal elements.